Role of network connectivity in intercellular calcium signaling 
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Abstract 

It is important to understand the coordinated performance of cells in tissue. One possible mechanism 
in this coordination involves intracellular Ca 2+ signaling. The topology of intercellular connections in 
tissue should also play an important role in this process. It is most relevant for plane tissues, in which 
the interaction between cells is due to gap junctions (epithelium, blood vessels). We demonstrate the 
importance of the topology of intercellular connectivity by investigating the properties of a model of 
Ca 2+ signaling for a small number of connected cells. 
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1 Introduction 

Calcium ions (Ca 2+ ) are one of the universal regulators of various intra- and intercellular 
processes jl|). In the presence of connections between cells, when each cell is a piece of tissue, 
Ca 2+ signal can be passed from cell to cell as intercellular Ca 2+ wave (J2I). The importance 
of intercellular Ca 2+ signaling (among other secondary messengers) in the coordinated per- 
formance of cells in tissues is currently under investigation. These waves have been seen in 
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a variety of cell types (0), (Q), (0), (Q). Intercellular Ca 2+ waves can be initiated by a focal 
mechanical, electrical, or hormonal stimulus. Often such waves travel from cell to cell in an 
oscillatory manner, thus forming periodic intercellular calcium waves (3), (0). 

Intercellular Ca 2+ signaling may occur through different pathways, for example through 
paracrine or junctional routes. In recent study (0) it was shown for astrocytes that prop- 
agation of intercellular Ca 2+ waves in glia is essentially determined by both gap junctions 
between cells and paracrine signaling by ATP, other nucleotides and/or their metabolites. 
A detailed theoretical and experimental study of calcium oscillations in astrocytes doublets 
and linearly coupled triplets is given by Ullah et al. (10). A model of intercellular calcium 
waves based on the action of calcium as a second messenger in paracrine signal transduction 
via calcium sensing receptor was proposed by Gracheva and Gunton (jlll ) and extended for 
a linear chain of cells by Kepseu and Woafo (1121 ). 

Liver, blood vessels, brain and epithelium are very good examples of tissues in which cells 
are connected with each other by gap junctions. It was shown by experiments and modeling 
that mechanically stimulated intercellular Ca 2+ waves in epithelial, glial and endothelial 
cells can result from the diffusion of inositol 1,4,5-trisphosphate (IP3) through gap junctions 



13l ). (1 14 ). (Il5l ). Although Ca + can move through gap junctions, a cellular Ca + wave 
is the cellular response to a diffusing of IP 3 , passing between cells via gap junctions and 
consequently participating in the release of Ca 2+ ions from the intracellular stores. However, 
in contrast with these cell types, the amount of IP3 diffusing through the gap junctions in 



connected hepatocytes is insufficient to induce a Ca + response 



(0), (0), 



il7|), yj). Thus, 



models based on the passive diffusion of IP3 between adjacent cells through gap junctions 
(Il5l ) cannot account for apparent intercellular Ca 2+ waves in hepatocytes. Hoefer (Il9l ) has 
proposed an alternative mechanism of intercellular calcium waves in hepatocytes, based on 
the diffusion of cytosolic calcium through gap junctions. Stochastic versions of both the (1171 ) 
and (19) models are given by Gracheva et al. (0). Clair et al. (16) have also shown that 
an agonist receptor gradient could be responsible for the apparent undirectional Ca 2+ waves 
observed in multiplets or in perfused intact livers. In this case, the main factor influencing 
Ca 2+ wave propagation in hepatocytes is the cellular distribution of agonist receptors. 

In addition, Tsaneva-Atanasova et al. (12 ll ) have shown for pancreatic acinar cells that 
intercellular calcium diffusion is necessary and sufficient to synchronize Ca 2+ oscillations 
in neighboring cells, whereas the function of intercellular IP 3 diffusion is to provide the 
propagation of oscillating Ca 2+ waves in tissue. 

In all the above mentioned models the junctional permeability was chosen to be constant. 
However, since calcium ions can close gap junctions (j22l ). (|23l ). this would appear to pose an 
obstacle to the propagation of intercellular Ca 2+ waves via gap junctions. Fortunately, this 
apparent paradox can be resolved by considering the temporal aspects of these responses. 
It is known that extended exposure to high concentrations of Ca 2+ (> 10 fiM) results in 
the closure of gap junctions (1221 ) . (1231 ) . It also appears that physiological concentrations 
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of Ca 2+ can reversibly close gap junctions but this response takes about 30 s (1241 ). Thus, 
because Ca 2+ waves propagate faster than the proposed closure times I25I ). the two possible 
intercellular messengers (Ca 2+ and IP3) (1261 ). (1271 ) can diffuse to adjacent cells to propagate 
the wave before Ca 2+ closes the gap junctions. Therefore, it would seem that a more accurate 
description of the gap junction permeability requires taking into account its dependence on 
the cytosolic calcium concentration, rather than treating it as a constant as in the earlier 
models (0), (0, (0, (Q- 

We believe that an important issue in modeling intercellular calcium signaling is to account 
for the cellular topology of the tissue. As a rule, flat tissues are modeled in the form of a 
two-dimensional square grid (1281 ). (I29l ). However, in reality, the structure of even flat tissues 
is significantly more complicated. Even small areas of tissue can exhibit different topological 
types of cell connections between each other (0), (0). In addition, in the case of local 
stimulation in which only one cell (not the entire tissue) is stimulated, it is important to 
know which cell is stimulated (for the same structure of cell connections). Thus, the aim of 
our study is to investigate the relation between the topology of cellular connectivity in flat 
tissues and the corresponding intercellular calcium signaling. 

Since gap junctions ensure a strict topology of intercellular connections, we take into ac- 
count in our study only tissues in which cells are connected with each other by gap junctions 
and the main mechanism of intercellular Ca 2+ wave propagation is a diffusion of a second 
messenger through these gap junctions. We do not exclude the possible influence of extra- 
cellular pathways on the intercellular calcium signaling. However, for the sake of clarity and 
simplicity, in this work we omit taking into account any second messenger diffusion through 
the extracellular matrix, which could be relevant, for example, in astrocytic networks (|9|). 
In our model, we take into account only the diffusion of IP3 through gap junctions as the 
main mechanism of Ca 2+ wave propagation in tissue (ll5l ). (12 ll ). (1251 ) . Also we include in 
our model the dependence of gap junctional conductance on concentration of cytosolic Ca 2+ 



2 The model description 



We use the minimal model (1301 ) as the basis of our model of intercellular Ca 2+ signaling. In 
this model (j30|) the total constant entry of Ca 2+ into the cytosol is uq + V\f3. This includes 
the influx uq from the extracellular medium and the IP3-stimulated Ca 2+ release v\(3 {f3 is 
the degree of saturation of the IP3 receptor). It is more convenient for us to include in the 
model an explicit dependence on the IP3 concentration in cell. Therefore we replace the 
IP 3 -stimulated Ca 2+ release V\j3 in ( j3ol ) by its linear dependence on IP 3 : ki n ip 3 IP^. 

We also include gap junction tunneling between nearby cells. In order to model Ca 2+ wave 
propagation the IP3 stimulus (which mimics the local application of a stimulus) is applied 
only to the first cell. We model this stimulus induced IP3 as following function of time for 
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the first cell: 

7 p 3 (— ).± <t<ti, 

jpistvmulus) = ^ IP $™*) it l< t < t 2j (1) 

where t - time, t\ - time moment by which IP 3 reaches its maximal value /p^ max ^ t 2 _ 
time moment at which IP3 stimulation is interrupted. Thus, we formulate the rate of IP3 
synthesis J syn t as follows: 

7 t, j-p(stimulus) /n\ 

'Jsynt ^rate 1 3 • ^^7 

Modeling the agonist generated IP3 concentration by such function is more convenient nu- 
merically than, for example, using a square impulse. In addition, we add to the model a 
term describing the degradation of IP3 (28) in each of the cells: 

_ VpIP^kp 



where k p is the half maximal rate of IP3 degradation (1281 ) . Without this term, the calcium 
levels reach unphysiological values in the asymptotic time limit (1251 ) . The values of all the 
parameters are give n in Table [TJ 

Lazrak et al. (124 ) have shown that during the fast increase of cytosolic Ca 2+ concentration 
in each of the coupled adjacent cells (from 0.2 fjM to 1 /xM), the gap junctional conductance 
decreases from 100 % to %. During the subsequent return of the level of cytosolic Ca 2+ to 
its normal value (about 0.2 fjM), the gap junctional conductance recovers with practically 
the same rate as its previous decrease (see Fig. 3 in (124)). This means that during each 
oscillation of cytosolic Ca 2+ in one of two adjacent cells, the gap junctions between these 
cells are transiently closed. However, this closure does not occur instantaneously, but during 
an interval of about 30 s. Relying on these data (Fig. 3 in (124) ) . we model the fast changes in 
the gap junctional conductance 7^ between adjacent cells i and j by the following function: 

k n - 

(4) 

Ua cytMax(ij) + K G 

where Ca cyt M ax (ij) is the larger of the two cytosolic calcium concentrations in cells i and j. 
Thus, 7y = 1 corresponds to a 100 % conductive gap junction between cell i and j. For 
other parameters see Table [TJ The gap junctional conductance in Fig. 3 in (0) decreases 
very rapidly, even with a small increase of cytosolic Ca 2+ . Therefore we choose n = 5 in the 
function (J3]) (see Fig. 1). The diffusion of IP3 provides the mechanism of intercellular Ca 2+ 
wave propagation in our model. 

Thus, the evolution of our system is described by the following differential equations (cf. 
10) , for parameter values see Table [TJ) : 

dC Clcyt(i) 

^ = kin + ki n jp 3 I i"3(i) — k ou tC 1 d C yt(i)~\~ 
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+ (k ,Ca u-, + k , Ca v°° l V CaA ^ \ h Ca2cyt{i) (*>) 

{ P °° l P0 ° ZW "Kl + Cal ool(i) K\ + Ca% t{i) } K* + Ca\ m ' (5) 



dCa voo i(i) _ ^ a %t(n 

K serca 



dt — Ca K! + Cal ym 

, LjU pool(i) oa cyf(») _ , „ , s 

rel Kl + Cal oom Kl + Ca% m W^M> W 

= Jsynt — Jdeg + k gap %j(IP^j) — IP3(i)) + k gap % k (I P 3 (^ k ) — IP3(i)), (7) 

where the indices i, j and k denote the cell number, and k gap is the coefficient of gap 
junctional communication for IP3 (see Tabled]). 

Thus, the index triplets (ijk) = (120) = (210) correspond to just two coupled cells with 
only ki2 = k 2 \ ^ 0. The case in which (ijk) = (100), Ca cyt Max(ij) = Ca cyt (i) and all 
junctional coefficients are equal to zero corresponds to just one cell. This is equivalent to 
equations for intracellular Ca 2+ signaling. It is also worth noting that in the case of index 
pairs (ijk) = (123) = (213) = (312) there are three possible different configurations of 
three coupled cells. If A43 = 0, we have a one-dimensional chain of cells, with the first cell 
interacting only with the second cell, the second cell interacting with both the first and 
third cells and the third cell interacting only with the second cell (see Fig. 2 A). The same 
linear chain with other orderings of the cells is shown in Fig. 2 B. If none of the junctional 
coefficients vanish, each cell interacts with two other cells (see Fig. 2 C). The right side of 
Fig. 2 shows graphs that correspond to these three configurations. 

The system can be easily generalized to the case of more than three cells. For example, all 
the configurations in two dimensions of four coupled cells are shown in Fig. 3. All these are 
biologically realistic in two dimensions. Examples of real cellular configurations in vitro are 
shown in Fig. 1 and 2 in (27) and Fig. 9 and 11 in (Q). The system of cells 1, 2, 3, 7 in Fig. la 
in (1271 ) (if the rest of the cells are removed) corresponds to the topological configuration in 
Fig. 3 A (cell 3 or 7 should be stimulated) and 3 B (cell 1 or 2 should be stimulated) in our 
study. The system of cells 1, 2, 5, 7 in Fig. la in (1271 ) (if the rest of the cells are removed) 
corresponds to the topological configuration in Fig. 3 C (cell 5 should be stimulated), 3 D 
(cell 2 or 7 should be stimulated) and 3 E (cell 1 should be stimulated) in our study. The 
system of cells 1, 3, 5, 7 in Fig. la in (1271 ) (if the rest of the cells are removed) corresponds to 
the topological configuration in Fig. 3 F (cell 1 should be stimulated) and 3 G (cell 3, 5 or 7 
should be stimulated) and the system of cells in Fig. 9 in (8) corresponds to the topological 
configuration in Fig. 3 F (cell B should be stimulated) and 3 G (cell A, C or D should be 
stimulated) in our study. The systems of cells in Fig. lc and Fig. 2 in (27) and in Fig. 11 
in (0) correspond to the topological configuration in Fig. 3 I and 3 J in the present paper. 
Finally, we believe that it is possible to extract the configuration shown in Fig. 3 H from a 
variety of cellular tissues. 
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3 Results 



We first consider the case of two cells. The first cell is stimulated from s to = 14 s 
with ip( max ^ — 4 2 /j,M applied at t\ = 4.2 s (see Fig. 4). Some IP3 from this impulse 
dissociates as Jd eg (eq. (J2D)- The remainder of IP 3 diffuses into the second cell through gap 
junctions. The concentration of IP3, which is received as the result of all these processes in 
the first cell, is shown on Fig. 5 (IP^m from s to 14 s). The concentration of IP3 which 
diffused into the second cell as the result of diffusion through gap junctions from the first 
cell is shown in Fig. 5 (IP^ 2 ) from s to 14 s). The change of IP3 in both cells occurs in 
this period by large advances, rather than smoothly. This is connected with the fact that 
the gap junctional conductance between cells is changed according to 7^ (eq. (j3J), rather 
than remaining constant all the time. Thus, as gap junctions are completely closed at a 
high concentration of cytosolic Ca 2+ in any two cells, the diffusion of IP3 through them 
temporarily ceases. 

Cytosolic Ca 2+ oscillations in the model cell occur at concentrations of IP3 about 1 /iM. 
This concentration of IP3 is the required level in the second cell that occurs at a time later 
than in the first cell. Therefore, the cytosolic Ca 2+ oscillations in the second cell arise after 
a certain time delay (see Fig. 6). We find in our model an inverse dependence of this time 
delay on the amplitude of the impulse of IP 3 which is applied to the first cell (see Fig. 6). 

For more than two connected cells, our model shows that the propagation of calcium waves 
depends on the spatial configuration of cells; this may be interpreted as the first indication of 
heterogeneity in intercellular Ca 2+ wave propagation. There are three possible configurations 
of interactions between three cells (see Fig. 2). If the location of the two neighboring cells is 
symmetric with respect to the first cell (Fig. 2 B, C), than one can see similar oscillations 
in both the second and third cells, following an initial delay of about 2 s from the onset of 
oscillations in the first cell. In the case of a linear chain (Fig. 2 A), there is a difference in 
the calcium oscillations in the second and third cells. There is not only a delay time between 
the onset of calcium oscillations in the first and second cells, but also between the second 
and third cells (both dela y ti mes are about 6 s). The value of this delay time slightly differs 
from that found in (|25l ). (1281 ) . This is due to the different choice of parameter values in the 
two studies. 

We also briefly describe the results for the case of four connected cells. The results for 
all configurations of four cells are summarized in Tables 2-4. We note that we study 
only in-plane (two-dimensional) cell configurations. It is possible to distinguish five different 
configurations of four cells on a plane. However, for the same topological configuration 
of cells, the Ca 2+ signaling greatly depends on which cell in the structure is stimulated. 
Taking into account this fact, we identify ten cell configurations. As in the case of three 
cells, the dynamics of Ca 2+ in cells symmetrically positioned with respect to the stimulated 
cell is identical; for example, see configurations "A", "B", "C", "E", "F", "G" and "H" in 
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Fig. 3. This holds when all the cells in these configurations are indistinguishable. Also, this 
assumption allows us to study the dependence of the intercellular Ca 2+ signaling only on the 
topology of the cell to cell connections in a given cell configuration. 

We now discuss the key features of signaling in four cells configurations. For example, in 
the case of configuration "H" (Fig. 3 H) the first cell has two connections with adjacent cells 
(there are four connections between cells altogether) and we find very good calcium wave 
propagation between cells (see Fig. 7). Also, the dynamics of Ca 2+ in cells 2 and 3, which 
are positioned symmetrically with respect to the stimulated cell 1, is the same (see Tables 2 
and 4, Fig. 7). The delay time between the onset of Ca 2+ oscillations in the stimulated 
cell and the beginning of oscillations in all other cells in the configuration depends on the 
ordering number of neighbors of a given cell. In the case of a direct connection with the 
first simulated cell (first order neighbors) the delay time varies from 2 s to 4 s (see Table 
3). For example, in configuration "I" cell 2 is the first order neighbor of cell 1. For a second 
order neighbor, there is one additional cell between the stimulated cell and a given cell. For 
example, in configuration "I" cell 3 is the second order neighbor for cell 1. For second order 
neighbors, the time delay for onset of Ca 2+ oscillation varies from 4 s to 11 s; this is larger 
then in the case of the first order neighbors (see Table 3). In addition, for cells next to 
the stimulated cell there is a difference in response depending on the presence or absence of 
connections with other second order neighbors (compare cell 3 and cell 4 in configurations 
"C" and "G", Tables 2-4). This suggests a heterogeneity in the Ca 2+ response in these 
cell networks. Of course, larger cell networks would have to be studied to strengthen this 
conclusion. For in-plane configurations of four cells there is only one neighbor of third order 
which is present in configuration "I". For the given parameters there are no oscillations in 
cell 4 of configuration "I", as indicated in Tables 2-4 and Fig. 8. 

Thus, each cell configuration can be characterized by a definite set of neighbors of each 
order; moreover, the propagation of intercellular Ca 2+ waves for each configuration depends 
on this set. According to the set of neighbors, all cell configurations can be divided into four 
classes: 

• one neighbor of first order, one neighbor of second order, one neighbor of third order 
(configuration "I"); 

• one neighbor of first order, two neighbors of second order (configurations "C", "G"); 

• two neighbors of first order, one neighbor of second order (configurations "A", "D", 
"H", "J"); 

• three neighbors of first order (configurations "B", "E", "F"). 

A fast decline of intercellular Ca 2+ waves is characteristic of the first two classes, whereas 
a steady expansion of intercellular Ca 2+ waves is characteristic of the last two classes (see 
Tables 2-4). For example, in the case of configuration "C" (Fig. 3 C) the first cell has only 
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one connection with another cell (there are four connections between all cells, as in "H") 
and this results in a decreased range of the calcium wave propagation from the first cell (see 
Fig. 7). There is an even stronger decrease in the range of wave propagation in the case of 
configuration "I", which is a linear chain of four cells (see Fig. 8). 

We now consider some interesting configurations of four cells in more detail. Configuration 
"C" (Fig. 3 C) belongs to the second class of configurations. In this configuration the cy- 
tosolic Ca 2+ displays the typical dynamics of this class, i.e., prolonged high-frequency Ca cyt 
oscillations in the first cell and shorter Ca cy t oscillations with lower frequency in the other 
cells. Cells 3 and 4 are equally spaced with respect to the first cell in configuration "C"; 
therefore their cytosolic Ca 2+ dynamics is identical. At the same time the total intensity 
of oscillations decreases with the distance from the first cell. Configuration "D" (Fig. 3 D) 
belongs to the third class in our classification. The topological difference between configura- 
tions "D" and "J" consists in the presence (configuration "D") and absence (configuration 
"J") of a connection between cells 2 and 3. This difference is manifested in Ca 2+ dynamics, 
primarily through small temporary distinctions (see Tables 3 and 4) and also in a weaker 
response in cell 2 in configuration "D", as compared with configuration "J" (see Table 4). 
This is connected with the fact that cell 2 in configuration "D" exchanges IP3 with cell 3 
in addition to cell 1, whereas there is no such exchange in configuration "J". Note that 
although in configuration "D", the disposition of cells 2 and 3 in the absence of cell 4 is 
symmetrical relative to cell 1 (Fig. 3 D), the cytosolic Ca 2+ dynamics in cells 2 and 3 is 
different (see Fig. 7). This is the result of the fact that cell 4 is an additional neighbor of 
cell 3, and an additional exchange of IP3 between these two cells takes place. That exchange 
decreases the duration of oscillations in cell 3 in comparison with cell 2 in configuration " D" 
(see Fig. 7). 

We find one more interesting regularity in the case of different configurations for four 
cells. We calculate the sum of the durations of Ca 2+ oscillations (last column in Table 2) 
and the sum of the number of peaks of Ca 2+ oscillations (last column in Table 4) in all 
four cells, for each configuration. Then we divide this total number of Ca 2+ oscillations by 
the corresponding total duration of Ca 2+ oscillations for each configuration and determine 
the mean frequency of oscillation for each configuration (Fig. 9). There is obviously a good 
correlation of this mean frequency with the classes of configurations described above. In 
other words, we find a correlation of each configuration class with a mean frequency for this 
class. This regularity allows us to conclude that the topology of the graph of intercellular 
connections has a more important influence on intercellular Ca 2+ signaling than one might 
think at first glance. 

We find that a variable gap junctional conductance can significantly influence the range 
of intercellular Ca 2+ wave propagation. This is most clear for the chain of four cells (con- 
figuration "I", Fig. 3 I). If the gap junctional conductance is described by 7^ (eq. (j4])), 
Ca 2+ oscillations take place only in cells 1, 2 and 3, but the Ca 2+ wave does not reach cell 4 
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(Fig. 8). If the gap junctional conductance is constant and independent of Ca cyt (7^ = 1), 
Ca 2+ oscillations take place in all four cells of configuration "I" (Fig. 10). Therefore, even 
a brief closure of the gap junctions between cells has a strong effect on the propagation of 
intercellular Ca 2+ waves. 

We also consider the effect of varying the parameter ka on the range of intercellular Ca 2+ 
wave propagation. If one increases ka, the gap junctional conductance 7^ decreases. For 
configuration "I" with ka = 0.45 fj,M, such a decrease of 7^ results in an abrupt limitation 
of intercellular Ca 2+ wave propagation. In this case, more intensive Ca 2+ oscillations with 
increased frequency and duration in comparison with ka = 0.4 fiM (see Fig. 8) take place 
in the first cell. In the second cell the frequency of Ca 2+ oscillations became appreciably 
smaller in comparison with ka = 0.4 fiM, whereas in the third and fourth cells there is 
only a negligible increase in cytosolic Ca 2+ above its baseline value. If one decreases ka, 
the gap junctional conductance 7^ increases. For configuration "I" with kc = 0.35 fjM the 
corresponding increase in 7^ results in an increase in the range of Ca 2+ waves and higher 
frequency oscillations similar to 7^ = 1 at ka — 0.4 fj,M (Fig. 10). 

Finally, we find the following results for all the configurations. First (denoted as A be- 
low), with a constant gap junctional conductance (7^ = 1), the number of Ca 2+ oscillations 
and their total duration increase, in contrast to the case with a variable ~fij(Ca cyt ). Second 
(denoted as B below), several things remain the same whether or not the gap junction con- 
ductance varies with concentration. These include the response of symmetrically positioned 
cells, the general dependence of frequency and duration of Ca 2+ oscillations in a cell on the 
number of its direct neighbors, and the relation between the time delays from the onset of 
oscillations in the first cell to the first oscillations in the other cells of the configuration. 

The above-mentioned features are the consequence of IP 3 diffusion through the gap junc- 
tions. At the same time a variable gap junctional conductance influences only the absolute 
(rather than relative) values of IP 3 . Thus, independently of whether there is a variable 
or constant gap junctional conductance, the qualitative dynamics of cytosolic Ca 2+ during 
intercellular signaling for different topological configurations of cellular connections remain 
the same. This conclusion remains true if a different model of intracellular Ca 2+ is used. 
Thus, the topology of cellular connections determines the signaling. 

The first features (A) described above are the consequence of the different absolute values 
of IP 3 which are acquired by each cell in configurations with open and closed gap junctions. 
In the case of the features (B), it is not important how much IP3 has arrived in a cell, since 
these are relative characteristics, and only the relative amount of IP3 which moves between 
neighboring cells is important. For example, it is important whether the neighboring cells 
obtained equal amounts of IP 3 or not; if not, the difference in the amounts is important. 
Thus features (A) are absolute characteristics of the oscillations, and are determined by the 
behavior of the gap junctions, whereas features (B) are relative characteristics of the values 
in cells (of a given configuration) with respect to each other and are determined by the 
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topological structure. 

Next, we plot in Fig. 11, the same diagram as in Fig. 9, but for the case of 7^ = 1. We 
find that we can distinguish only the first ("I") and the second ("C", "G") configuration 
classes, but we cannot distinguish between the third and fourth classes in our configuration 
classification. To conclude, a variable gap junctional conductance has a modulating influence 
on intercellular Ca 2+ signaling by reinforcing the role of the cellular connection topology, 
but leaves the qualitative aspects unchanged. 

We have also considered Ca 2+ diffusion through the gap junctions between cells, in addition 
to IP3 diffusion between cells. We thus added an exchange term to eq. §5§ in the form 
k ga pjij(Ca cyt Q) —Ca cyt (i)), similar to eq. ([7j). The results we obtain allow us to conclude that 
there are two separate mechanisms involved in intercellular signaling. First, IP3 diffusion 
provides the mechanism for calcium wave propagation between cells, as was proposed in 



(I28l ). Second, the diffusion of Ca 2+ , is necessary (in addition to IP3 diffusion), to synchronize 
the oscillations among the neighboring cells. Similar conclusions were drawn by Tsaneva- 
Atanasova et al. (0), which are in complete agreement with this work. 



4 Conclusion 

Intercellular Ca 2+ signaling is a complex and nonlinear process; its mechanism depends on 
the kind of tissue and the method used to stimulate the cells. In this work we study only 
plane tissues, cells of which are connected to each other by gap junctions, such as epithelium 
and blood vessels. We suggest that the diffusion of IP3 through gap junctions is the main 
mechanism for intercellular Ca 2+ wave propagation in such tissues. That is, we consider 
intercellular topological structures which are common to epithelium and endothelial tissues, 
but we do not model the details of the intracellular signaling specific to these cell types. 
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cell cultures arises as the 



The intercellular Ca 2+ wave propagation in airway epithe 
result of mechanical stimulation of a single cell in the tissue (1251 ) . The front of this wave is 
heterogeneous in different directions; thus one can observe many kinds of "branches", "paws" 
and so on, as well as an oblongness of the front on a whole in some directions (1251 ) . The 
reasons for such a heterogeneous wave propagation might include both the heterogeneity of 
the tissue structure (for example, the different number of gap junctions in the boundaries 
between adjacent cells, the different number of neighbors for each cell, the different shape 
of cells and so on), and the internal heterogeneity of the cellular structure (I.V.Dokukina, 
M.E.Gracheva, E.A.Grachev - unpublished results). In the current study we show that differ- 
ent numbers of neighbors results in a heterogeneity of intercellular Ca 2+ wave propagation, 
while varying the gap junctional conductance amplifies this influence. 

We believe that the topology of intercellular connections in tissue has an essential influence 
on Ca 2+ signaling. We demonstrate this influence for the case of a small (three of four) 



number of cells, using the minimal model (1301 ) for the description of intracellular processes 
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In addition, we take into account the dependence of the gap junctional conductance on 
the cellular concentration of cytosolic Ca 2+ . Even in the simplest case of three cells we 
find differences in Ca 2+ signaling between the linear chain of three cells (Fig. 2 A) and 
the closed chain in which each cell is connected to its adjacent cells (Fig. 2 C). We find a 
stronger influence of the topological structure of intercellular connections on Ca 2+ signaling 
for different configurations of four cells. We classify all possible plane configurations of four 
cells in four different categories, depending on the number of neighbors of a definite order. 
We find that the type of intercellular Ca 2+ signaling in each configuration depends on the 
class of this configuration. Furthermore we find that by knowing the mean frequency of Ca 2+ 
oscillations in any configuration, it is possible to refer this configuration to one of the four 
categories. 

Our results suggest that during any investigations of intercellular Ca 2+ signaling in a 
particular tissue, it is necessary to pay attention to the possible influence of the topology 
of cellular connections. This is most important for small pieces of tissue, because in this 
case the influence of topology is most evident. Note that we consider only two-dimensional 
configurations of cells in this paper. However, it would be interesting to consider the role of 
topological structure in three-dimensional tissues, which is a subject of our future research. 



Also, Tsaneva-Atanasova et al. (12 ll ) have shown that an analysis of the point-oscillator model 
is insufficient to determine which kinds of oscillations will appear in clusters of pancreatic 
acinar cells with realistic geometries. They also find that geometrical factors play a crucial 
role in both intra- and intercellular Ca 2+ signaling. Therefore, it is important to carry out 
a study of spatially distributed models of the most interesting topological configurations of 
cellular connections. 

While modeling the intercellular calcium signaling it is necessary to recognize that the 
transition from one cell to tissue is not a simple quantitative conversion. It is not a transition 
from one cell to a chain of cells or some grid of cells; rather, it is a more complex transition 
defined by the topology of tissue. We have shown that the dynamics of tissue in small areas 
(and also in larger regions - I.V. Dokukina, A. A. Tsukanov, M.E. Gracheva, E.A. Grachev, 
unpublished results) is defined not by a chain or grid of cells, but by the topology of the 
connection graph of cells with each other. It is possible to judge the tissue topology from 
the dynamics of signaling; in turn, the tissue topology can change the intercellular signaling 
dynamics, enforcing its own rhythm. Thus, the intercellular dynamics of calcium and the 
corresponding topology of cell connections form an unbroken, complete system. 
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Figure Captions 

Fig. 1. Behavior of function 7^ = ^-^ — vpr for different values of the power n. Here 

cytM ax(ij) G 

kij = 1 /iM at n — 5, fcjj = 1.84 /xM at n = 3 and fcy = 3.94 /xM at n = 2. We represent 
7y in percentile form, or in other words normalize 7^ by max (7^-) to obtain the percentile 
value of 7y from its maximum. For other parameters see Table HJ 

Fig. 2. Different structures of connections between three cells. The grey color of cell 
1 means that this cell is stimulated. The right side shows graphs that correspond to each 
configuration. 

Fig. 3. Different structures of connections between four cells. The grey color of cell 1 
means that this cell is stimulated. The right side shows graphs that correspond to each cell 
configuration. 

Fig. 4. Schematic representation of stimulus induced IP3 as a function of time. 

Fig. 5. Intercellular calcium signaling between two cells based on diffusion of IP3 through 
gap junctions. Gap junctional conductance depends on cytosolic Ca 2+ . For all parameters 
see Table [U 

Fig. 6. The dependence of the time delay between initiation of oscillations in the first 
(stimulated) and the second cells on the amplitude of the initial impulse of IP 3 in the first 
cell. 

Fig. 7. Intercellular calcium signaling between four cells of configurations "C", "D" and 
"H" (Fig. 3 C, D and H) based on diffusion of IP3 through gap junctions. Gap junctional 
conductance depends on cytosolic Ca 2+ . For all parameters see Table [Q 

Fig. 8. Intercellular calcium signaling between four cells of configuration "I" (Fig. 3 I) 
based on diffusion of IP3 through gap junctions. Gap junctional conductance depends on 
cytosolic Ca 2+ . For all parameters see Table HJ 

Fig. 9. Axis X: all configurations for four cells are separated into four different categories, 
taking into account the number of neighbors (see text for details). Axis Y: the mean fre- 
quency of oscillation for all cells of each configuration (mean frequency = ^ Z^tZ of a11 oscillations) 
with ^ij{Ca cy t). 

Fig. 10. Intercellular calcium signaling between four cells of configuration "I" (Fig. 3 
I) based on diffusion of IP3 through gap junctions. Gap junctional conductance does not 
depend on cytosolic Ca 2+ . For all parameters see Table [TJ 

Fig. 11. Axis X: all configurations for four cells are separated into four different categories, 
taking into account the number of neighbors (see text for details). Axis Y: the mean fre- 
quency of oscillation for all cells of each configuration (mean frequency = duration °^ an oscillations) 
with 7^ = 1. 
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Table 1: Model parameters 



Designation 


Parameter values 


kin 


1 (jM ■ s- 1 


kout 


6 s- 1 


kin IP 3 


1 s- 1 


kpool 


Is" 1 


k re i 


500 fM • s" 1 


kserca 


65 (M ■ s- 1 


K x 


2 nM 


K 2 


0.9 fjM 




1 fM 


v P 


0.01 s- 1 


hp 


1 


kij 


1 (M 


k G 


0.4 (M 


kgap 


0.005 s- 1 


krate 


Is- 1 
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Table 2: Duration of Ca 2+ oscillations (s) 



Configuration 


Cell 1 


Cell 2 


Cell 3 


Cell 4 


Total 


A 


19 


14.5 


14.5 


15.8 


63.8 


B 


18.7 


14.6 


14.6 


14.6 


62.5 


C 


19.7 


15 


10.3 


10.3 


55.3 


D 


17.6 


16.7 


14.6 


13.6 


62.5 


E 


18.7 


14.6 


14.6 


14.6 


62.5 


F 


18.7 


14.6 


14.6 


14.6 


62.5 


G 


19.7 


15 


10.3 


10.3 


55.3 


H 


19 


14.5 


14.5 


15.8 


63.8 


I 


21.4 


15.4 


7.5 




44.3 


J 


16.9 


18 


16.4 


11.4 


62.7 
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Table 3: Time delay from the beginning of Ca 2+ oscillations in the first cell (s) 



Configuration 


Cell 2 


Cell 3 


Cell 4 


A 


2.5 


2.5 


4 


B 


2 


2 


2 


C 


4 


8.5 


8.5 


D 


2 


2.5 


5.5 


E 


2 


2 


2 


F 


2 


2 


2 


G 


4 


8.5 


8.5 


H 


2.5 


2.5 


4 


I 


3 


11 




J 


2 


3 


7.5 
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Table 4: The number of Ca 2+ oscillations (items) 



Configuration 


Cell 1 


Cell 2 


Cell 3 


Cell 4 


Total 


A 


9 


6 


6 


6 


27 


B 


8 


6 


6 


6 


26 


C 


13 


6 


4 


4 


27 


D 


9 


7 


6 


5 


27 


E 


8 


6 


6 


6 


26 


F 


8 


6 


6 


6 


26 


G 


13 


6 


4 


4 


27 


H 


9 


6 


6 


6 


27 


I 


15 


7 


3 




25 


J 


9 


8 


6 


4 


27 



19 



100 
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